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ABSTRACT 

We develop statistical methods for identifying star streams in the halo of the Milky Way galaxy 
that exploit observed spatial and radial velocity distributions. Within a great circle, departures of 
the observed spatial distribution from random provide a measure of the likelihood of a potential star 
stream. Comparisons between the radial velocity distribution within a great circle and the radial 
velocity distribution of the entire sample also measure the statistical significance of potential streams. 
The radial velocities enable construction of a more powerful joint statistical test for identifying star 
streams in the Milky Way halo. Applying our method to halo stars in the Hypervelocity Star (HVS) 
survey, we detect the Sagittarius stream at high significance. Great circle counts and comparisons 
with theoretical models suggest that the Sagittarius stream comprises 10% to 17% of the halo stars 
in the HVS sample. The population of blue stragglers and blue horizontal branch stars varies along 
the stream and is a potential probe of the distribution of stellar populations in the Sagittarius dwarf 
galaxy prior to disruption. 

Subject headings: galaxy: kinematics and dynamics — Galaxy: structure — Galaxy: halo — Galaxy: 
stellar content — galaxy: individual (Sagittarius dwarf galaxy) — stars: horizontal- 
branch — stars: blue stragglers 



1. INTRODUCTION 

The timescale for halo stars to exchange their energy 
and angular momentum is long compared to the age of 
the Milky Way galaxy. Thus, the stellar halo of the Milky 
Way provides a unique laboratory for studying hierar- 
chical galaxy formation. Theoretical simulations show 
that the remnant debris of hierarchical galaxy formation 
should be visible as st ar streams in the Mi l ky Way halo 
(I Johnston et al.lll996t iHardi ng et a" Lll2Q01l: [Abadi et al.l 
[2003; Bul lock fc JohnstonI 12005: Fo nt et al.l I2OO60 . In- 
deed, observers have discovere d the tidal de bris of the 
Sagittarius dwarf gala xy (Ibata et al .l 11994*) wrapping 
around the Milky Way (jMajewski et al. 2003) and other 
stellar overdensities t hroughout the stellar halo (e.g., 
iBelokurov et al.l [20061 ). The vast majority of these de- 
tections come from star counts without radial velocity 
measurements. 

We describe a powerful technique for identifying tidal 
star streams that exploits the additional information 
provided by a complete radial velocity survey. With 
large stellar radial v elocity surveys suc h as RAVE 
'S tcinmctz et al. 2001 iZwitter et a"l1l2008D and SEGUE 



Yannv et al.l 120091 now available, approaches including 
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radial velocities enable better constraints on structure in 
the halo. 

To identify statistically significant tidal structures, we 
examine the surface number distribution and the ra- 
dial velocity distribution of stars in great circles. Both 
the surface density and the radial velocity distributions 
probe the origin of structure in the halo. A radial velocity 
test for structure applies independently of the underlying 
spatial distribution. Our technique opens up the possi- 
bility of finding tidal star streams based on structure in 
velocity space. Combining our spatial and radial veloc- 
ity statistical tests provides a more powerful joint test 
for identifying star streams. 

Previous studies of positions and velocities have 
focused o n known tidal st r eams, such as Sagittar- 
ius (e.g., lYanny et all 120091: iRuhland etall 1201 Ih . or 
have analyzed disconnected p a tches scattered across th e 
sky (" Starkenburg et al.l 120091: iSchlaufman et al.l I2009D . 
Known tidal debris, however, is contiguo us over areas of 
hundreds or thousands of square degrees ([Johnston et al.l 
■ 1999; Jbata et al., 2001) . 

As a first demonstration of our technique, we search for 
structure in the HVS radial velocity survey, a complete, 
non-kinematically selected sample of late B-type stars 
covering a contiguous area of more than 8400 deg^ on the 




Fig. 1. — Solid angle, f2, subtended by a great circle segment 
with arc length 9, width A, and pole at galactic coordinates (Z, b). 



sky (IBrown et alJ [20051 12nn6bl . WmA [20Toh . The HVS 

survey targets blue stragglers and blue horizontal branch 
(BHB) stars within the halo and is distinct in its uniform 
velocity data for a well-defined sample of potential halo 
stars. 

We detect significant structure in the sample of halo 
stars within the HVS survey. Most of this structure is 
attributed to the Sagittarius stream. Wc compare our de- 
tections with the Law & Majewski (2010) A^-body model 
for the Sagittarius dwarf, derive estimates for the fraction 
of HVS stars within the Sagittarius stream, and assess 
the relative fractions of blue stragglers and BHB stars 
along the stream. 

In 321 we describe the statistical technique. In ^ 
we apply this technique to the HVS survey. In §4l 
w e compare the survey st ars with the A^-body model 
of iLaw fc Maiewskil ()2010[ ) in order to explore the na- 
ture and surface density of stars within the Sagittarius 
stream. We conclude in fJSl 

2. TECHNIQUE FOR IDENTIFYING STAR STREAMS 

Identifying structures in the Galactic halo is a challeng- 
ing statistical problem. The structures are low contrast 
and may extend across the entire sky. For most stars 
in the Galactic halo, we have ready access to only half 
of phase space. Radial velocity and angular position are 
the only robust coordinates available f or large samples of 
stars more distant than ^ 10 kpc (e.o.. lSchlaufman et al.l 
l2009l) . 

Methods of identifying star streams and other sub- 
structures in the Milky Way halo include pole counts 
(iLvndcn-Bell & Lvnden^Bell IMl) , great circle counts 
(Johnston et al. 1996), star count maps (e.g.. llbata et al.l 
[2002; Belokurov ct al. 2006^ group fi nding algorithms 
(jSharma et al., .2011; .Xue et al.l 120 111) , and combined 
analyses of the spatial a nd radial velocity distri butions 
of halo populations (e.g.. iSchlaufman et al.|[2009[ ). Most 
of these methods include an assessment of the statistical 



significance of the structures they detect. 

The method we develop contains elements of previous 
approaches refined for application to radial velocity sam- 
ples that are complete over some substantial portion of 
the sky. Our approach has three steps: (1) we analyze 
great circle counts based on a nearly uniform, dense grid 
of sampling poles, (2) we compare the great circle radial 
velocity distribution with the radial velocity distribution 
for the entire survey, and (3) we combine the great circle 
counts and velocity distribution comparison into a single 
statistic. We estimate the significance level, or probabil- 
ity of false positives, among the structures we detect. 

2.1. Great Circle Counts 

The orbits of tidal debris from disrupted satellite galax- 
ies in the Milky Way halo persist for several Gyr. Be- 
cause the Galactic halo is approximately spherical, a 
satellite's orbit in the sky should approximate a great 
circle. The "fossil" debris shoul d be distributed along 
the orbit of th e satellite {e.g., [Johnston et al.l 119961 : 
[Bullock fc Johns ton 2005). 

[Johnston et al.l (11996) propose using star counts in 
great circles for detecting tidal debris against background 
halo stars. They define a great circle by its width and 
by the position of its pole relative to the Galactic pole. 
The pole is the unit normal to the plane of the great cir- 
cle. The latitude and longitude of the pole on the unit 
sphere then determine the location of it s corresponding 
great circle, as illustrated in Figured] [Johnston et al.l 
([1996( ) compare counts of stars observed within a great 
circle with the number of stars expected from a random 
distribution of stars on the sky. We generalize this ap- 
proach for an observational sample that covers only a 
portion of the sky. For simplicity, we refer to the num- 
ber of stars in a great circle segment as the great circle 
count. 

For a great circle segment of arc length 9 and width 
A within the sample region, the solid angle subtended 
the segment is ^1 = 29s\n{^) ~ OA for small A. The 
probability, p, that a star placed at random in the sam- 
ple region falls within the great circle segment is the ra- 
tio between the area of the segment, ft, and the area, 
^sampie^ of the entire sample region (Figure [1]) 
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for small A. (1) 



For n stars randomly distributed in the sample region, 
the number, k, of stars in the great circle segment follows 
a binomial distribution with expected value E{X) = np 
and variance Var{X) — np{l —p). The cumulative bino- 
mial distribution gives the probability, pgCj of observing 
k or more stars in the great circle segment 

PGC = Pr(X > fc) = Qp^" (1 - pf-' (2) 

j—k 

We evalu ate this expressio n using the incomplete beta 
function fiPress et al.l 11993) . This probability measures 
the statistical significance of each observed great circle 
count. 

We sample the survey region by placing poles on a 
lattice of Galactic longitude, I, and latitude, b, derived 
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from a nearly uniform distribution of antipodally sym- 
metric point s on the unit sphere using the method of 
iKoavl (|2011[ ). Note that a pole and its antipode rep- 
resent the same great circle. We then record observed 
counts for all segments of great circles corresponding to 
these poles that lie within the sample boundaries. 

To account for multiple statistical tests on the same 
set of data and to avoid overstating the statistical signif- 
icance of the findings, we m o dify t he great circle count 
technique of 'Johnston ct al.' (1996). Using the Bonfer- 
roni correction (Abdi 2007), wc adjust the probability of 
each individual test, pcc in equation to account for 
multiple tests by multiplying the probability by the total 
number of tests of the data. The Bonferroni correction 
protects against false positives but at the expense of more 
false negatives. The Bonferroni correction assumes that 
the tests are independent and represents a lower bound 
when they are not. Thus, our estimate of the overall 
probability of false positives is conservative. In fact, we 
expect that none of the highest confidence streams we 
detect in the HVS survey sample are false positives (see 

2.2. Radial Velocity Distribution 

Radial velocities of stars are a potentially important 
probe for the presence of star streams. The velocity dis- 
tribution for halo stars is nearly Gaussian with average 
line-of-sight velocity Vrf — km s~^ and velocity dis- 
persion CTjj_^^J^6 km_s^ ^ in the Galactocentric rest 
frame (jBrown et al.|[2010[ ). The radial velocity disper- 
sion for stars in a s t ream should be < 10 km s~^ (e.g., 
iHarding et all 120011: iMaiewski et~all I2004D . Relative to 
the large dispersion of the halo, this small velocity dis- 
persion should enable detection of streams in velocity 
space. 

We search for evidence of streams by comparing the 
distribution of radial velocities for the stars in a great cir- 
cle with the observed distribution of radial velocities for 
the entire sample. If the distribution of radial velocities 
on the sky is randomly drawn from this underlying distri- 
bution, the radial velocities of stars lying in an arbitrary 
great circle should follow the same distribution. Rather 
than assuming a Gaussian radial ve locity distribution of 
background halo stars a priori {cf. IHarding et aT1[200ll: 
ISchlaufman eFall [20091) . we use the observed radial ve- 
locity distribution for the sample. In this approach, the 
two-sample Kolmogorov-Smirnov test provides a non- 
parametric estimate of the probability, pv^f , that the ra- 
dial velocity distribution for stars within a great circle is 
randomly drawn from the same distribution as the HVS 
survey. We compute pv^^ for great circles with poles on 
the same grid we use for the great circle counts. For each 
segment of a great circle that lies within the boundaries 
of the HVS sample, we have pcc and pv^f ■ 

We combine the great circle count and radial velocity 
probabilities to form a joint test under the assumption 
that spatial position and radial velocity are independent. 
The probability of observing a group of stars and their 
radial velocities within a great circle segment of width A 
is then the product of the individual probabilities for each 
of these events, pjoint — Pgc x Pv^-f ■ Although position 
and radial velocity may not be truly independent, the 
combined probability gives a measure of the joint depar- 



ture of the spatial and velocity distributions from those 
for a set of points randomly distributed on the sphere 
and randomly drawn from the observed radial velocity 
distribution. 

2.3. Limitations of the Technique 

There are several important observational and theo- 
retical limitations to our ability to detect tidal streams. 
Observational limitations include the limited sky cover- 
age of a sample and the limited depth and density of 
a sample. Theoretical limitations include departure of 
the orbit of the stream from a great circle and mult iple 
orbits or wraps of the streams. IHarding et al.l ()2001l ) re- 
view many of the issues involved in detecting tidal debris 
in the Galactic halo. 

Existing radial velocity surveys sample only a fraction 
of the sky. Consequently, we can only access great circle 
segments with arc lengths 9 < 360°. The observed mag- 
nitude limit of a radial velocity survey defines the volume 
of the halo and the part of the orbit we sample. We are 
potentially less sensitive to the apocenter passage than 
to the pericenter passage of the orbit. More stars are 
stripped from the progenitor of a star stream near peri- 
center, probably enhancing our ability to detect streams 
and possibly leading us to overestimate the density of the 
stream. 

Our sensitivity depends on the width of the great cir- 
cle that we choose. The overall density of the sample 
sets a lower limit on the width of great circles we can 
profitably explore. As we increase the width of the great 
circles we use to sample the distribution, dilution of po- 
tential streams by background stars increases. This prob- 
lem is particularly serious for the intrinsically narrowest 
streams. 

For computational convenience we approximate orbits 
as great circles. However, the sun's offset of approxi- 
mately 8 kpc from the Galactic center means that orbits 
do not appear exactly as great circles on the celestial 
sphere, an effect that diminishes with the size of the orbit 
and the proximity of the sun to the plane of the orbit. 
Furthermore, the orbits depart from great circles for a 
nonspherical Galactic potential. 

A nonspherical Galactic pot ential causes precessio n 
of the orbit of the progenitor ( Law fc Ma iewski '2010^. 
Stars stripped during different passes of the progenitor 
may then appear in the same observed great circle. In 
our technique, one indication of this phenomenon is the 
appearance of multiple peaks in the observed radial ve- 
locity distribution for the great circle approximating the 
orbits. Multiple wraps may also increase the great cir- 
cle count. In addition, intersection with other streams 
or structures may produce multiple peaks in the radial 
velocity distribution as well as an increase in the great 
circle count. 

3. APPLICATION TO THE HYPERVELOCITY STAR 
SURVEY 

We demonstrate our technique by applying it to the 
spectrosc opic sampl e of distant halo stars from the HVS 
surv ev of iBrown et a l. (2005, 2006a,b, 2007a,b, 2009a,b). 
The lBrown et al.l (|2010 ) catalog of 910 late B-type stars 
is 93% complete over the magnitude range 17 < 50 < 20.5 
and covers more than 8400 deg^ of the Sloan Digital Sky 
Survey (SDSS) Data Release 6 (DR6) imaging footprint. 
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Fig. 2. — Galactic coordinates {I, b) of the 881 halo stars in the HVS program survey with |fr/| < 275 km s~^. As indicated in the legend 
to the right of the main figure, color represents radial velocity in the range —300 km s^^ to +300 km s~^. 



These halo stars are a mix of evolved blue stragglers and 
BHB stars. Because the selection is well understood, 
this survey and others that cover a large area of the sky 
uniformly are particularly powerful probes of the halo. 

3.1. Data 

iBrown et al.l (|201Clt ) measured radial velocities using 
the Blue Channel Spectrograph on the 6.5m MMT. They 
derived velocities by constructing appropriate stellar 
templates and then applying standard cross-correlation 
routines (Kurtz & Mink 1998). The typical radial veloc- 
ity uncertainty is ±12 km s~^. IBrown et al.l (|2010[ ) trans- 
form heliocentric velocities (vheiio) into Galactocentric 
rest frame velocities assuming a circular velocity of 220 
km s-i and a solar motion of {U, V, W) = (10.0, 5.2, 7.2) 
km s~^ (|Dehnen fc Binnevl [T998'). For each star, we use 
this Galactocentric rest frame velocity, Vrf. 

From the 910 stars in HVS survey, we construct a sam- 
ple of the 881 stars with \vrf \ < 275 km s~^. The upper 
limit on the velocity removes unbound HVSs. The 881 
stars in our sample range in heliocentric distance from 12 
to 91 kpc. Figure [2] shows the distribution of the sample 
stars on the sky; the color of the points indicates their 
radial velocities. 

3.2. Great Circle Counts 

For the HVS survey, we calculate pcc (equation ©) 
for great circles with width A = 5°. The surface den- 
sity of stars in the survey constrains A. We chose this 
width by experiment. The survey is too sparse to de- 
tect streams with A < 5°; adopting A > 10° dilutes 
the signal of streams. Other investigators use A = 5° 
( jMajewski ej al 2 003), 10° (I bata e t al. 2001) and 12° 
(|Ibata et al.ll2002f ). We distribute 2.4 x lO'* antipodally 
symmetric p oints nearly uniformly on the unit sphere 
(|Koavll20Tll) . approximately one pole per 1.7 square de- 
grees. Since the poles are antipodally symmetric, this 
represents 1.2 x 10'^ unique great circles. We consider 



only great circles that intersect the sample region with 
arc lengths 9 > 25°. For shorter arcs, small number 
statistics dominate the results. Our dense grid of sam- 
pling poles f §2.1|) produces 9884 such great circle seg- 
ments. Many of these segments overlap. 

Figures 3 and 4 show results for the great circle counts 
derived from the unadjusted pcc without the Bonfer- 
roni correction (Figure [3]) and the adjusted poc with the 
Bonferroni correction (Figure |4]). At the position of each 
pole, we encode log^Qjl/pGc} in the color scale indicated 
in the legend. We plot only results with p- values of 5% 
or less {pgc 0.05). The plot has an inherent symmetry 
as poles at antipodes, (/, b) and (/ -I- 180°, —b), represent 
the same great circle. 

Comparison of these two figures illustrates the effect 
of accounting for multiple statistical tests applied on the 
same data set. Applying the Bonferroni correction in- 
creases the probabilities that the observed great circle 
counts result from a random distribution of stars on the 
sky by a factor of 9884. The values of p'^^^*'^'^ Figure 

H range from 5.0 x 10"^ to 3.63 x 10"^. 

Using the adjusted pac, the three most significant 
poles in Figuregjlie at (103.1°, 12.4°), (104.5°, 12.4°) and 
(121.6°, 20.2°) with p^'*"'**'^'^ = 3.63 x 10"^ 4.99 x 10"^ 
and 1.10 X 10~^, respectively. Features of this signifi- 
cance have much less than a 1% chance of appearing in 
a random distribution. We list the 25 most significant 
poles and identify all of them with known structures in 
Table [T] The two most significant poles coincide with 
the Sagittarius stream (S|3]4] below). 

3.3. Using Radial Velocities 

Radial velocities provide additional information that 
aids in identifying star streams, in discriminating among 
them, and in distinguishing them from other structures. 
Figure [S] illustrates the power of comparing the velocity 
distribution for stars in a great circle with the corre- 
sponding distribution for the entire HVS survey. The 
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Fig. 3. — Significance, logio(l/pGC')i of the observed great circle 
count for pole at (I, b), where pec represents the unadjusted prob- 
ability that the observed count results from a random, isotropic 
distribution of stars, 6 > 25° and A = 5°. As indicated in the 
legend to the right of the main figure, color denotes significance in 
the range logio(l/PGc) = 1.5-12.5. 




Fig. 4. — As in Figure[3]for logj^Q(l/p' 
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green histogram compares the stellar radial velocity dis- 
tribution for the great circle with pole at (89.3°, 15.0°) 
with the velocity distribution for the complete survey 
(red histogram). For this, the most significant structure 
in radial velocity in Table[l] the two-sample Kolmogorov- 
Smirnov test yields a probability of py^^. = 1.47 x 10~^ 
that the radial velocities in the great circle were drawn 
from the same distribution as the entire sample. 

The radial velocity distribution within this great circle 
has two peaks, one at Vrf « —80 km and another at 
Vrf w —190 km s~^. These peaks do not appear in the 
distribution for the entire HVS survey sample. The ex- 
pected radial velocity dispersion within a stream should 
be < 10 km s~^; thus, these peaks suggest two distinct 
physical components in the velocity distribution. In H'dAl 
we identify this structure in velocity with the Sagitarrius 
stream. 

Using our HVS survey sample, we calculate the prob- 
ability, pv^f , that the radial velocities for stars in the 
entire sample and for stars within a great circle come 
from the same distribution for all great circle poles in 
our grid. Significant poles are located near the signifi- 
cant great circle count pole at (89°, 15°). 

The most significant pole based on radial velocity alone 
has pv^f = 8.31 X 10"^ appears at (186.6°, 2.0°), and 
corresponds to a near polar great circle. Figure [5] shows 
the velocity distribution of the stars in the great cir- 
cle corresponding to this pole. This detection probably 
results from two separate substructu res, the Sagitarrius 
stream and the Virgo overdensity (jNewberg et al... 20071) 
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□Great Circle {89.3°, 15.0°) (n = 1Q6) 



-220 -165 -110 -55 55 110 165 220 

Radial Velocity in Galactic Rest Frame, v (km s~^) 



Fig. 5. — Distribution of galactoccntric radial velocities, 



for 



sample stars in the great circle with pole at (Z, b) = (89.3°, 15.0°) 
intersecting Sagittarius stream branch A compared with those of 
entire sample. The two sample Kolmogorov-Smirnov test yields a 
probability Pv^f = 1-47 X 10~^ that these two distributions are 
drawn from the same distribution. 




-^5 -220 -165 -110 -55 55 110 165 220 275 

Radial Velocity in Galactic Rest Frame, (km s~^) 

Fig. 6. — Distribution of galactocentric radial velocities, v^j, for 
sample stars in the great circle with pole at {l,b) = (186.6°, 2.0°) 
intersecting the Sagittarius stream and the Virgo overdensity com- 
pared with those of entire sample. The two sample Kolmogorov- 
Smirnov test yields a probability Pv^f = 8-31 x 10~^ that these 
two distributions are drawn from the same distribution. 



as the great circle intersects both. Some stars may be- 
long to the Sagitarrius stream f H4.ip . Others may belong 
to the Virgo overdensity, which has two peaks in its ra- 
dial velocity distribution at Vrf = —49 km and at 
Vrf = —171 km (jVivas et al.ll2008l ) that match those 
of stars in this great circle. 

Next we combine both spatial and radial velocity in- 
formation in a joint statistic. Figure [7] shows P°fo^i'^'^ = 
Pgc"''*'^'* X PVrj the positions of all of the poles in our 

grid. We plot logio{l/Pjoint*^''}- Many features here are 
the same as in the map of significant great circle count 
poles (Figured]). 

We list the 25 most significant poles in Pj'^l^f'^'^ 
and their associations with the Sagitarrius stream and 
other known structure in Table [H The two most 
significant poles in Figure [7] and Table [T] correspond 
to branch A of the Sagittarius stream identified by 
iBelokurov et al.l ()2006[ ). These poles are also the most 
significant in p'^^"*^'^ alone and occur at (103.1°, 12.4°) 




Fig. 8. — Great circles in Galactic coordinates (l,b) corresponding to the 25 poles with most significant Pj^'i^J*'^'* = p^'^""*'^'* x Pv^f 
Solid gray lines mark great circles affiliated with the Sagittarius stream. Dashed gray lines denote great circles with superpositions of the 
Sagittarius stream and other structure. The large green dot labeled Sgr marks the location of the Sagittarius dwarf galaxy. As indicated 
in the legend to the right of the main figure, color denotes radial velocity in km . 
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Fig. 7. — As in Figure |3] for the adjusted joint probability, 

1 /-I I adjusted\ t 

log(l/Pjo,„t ), where p Joint ^ 



adjusted , , 
^ PGC ^ PVrf ■ 



adjusted 



4.63 X 10 



-10 



and 



and (104.5°, 12.4°) with pj' 

1.09 X 10~^, respectively. Although the radial velocity 
is helpful in selecting other Sagittarius poles, the values 
for pv^f alone for these two poles are of low significance 
(Table p. 

3.4. Detection of the Sagittarius Stream 

llbata et al.l ()2002D analyzed M giant source counts in 
26.4% of the sky using the 2MASS Second Incremental 
Data Release and identified a peak corresponding to a 
Sagittarius plane with pole at (95°, 13°) fo r A = 10° 
Using carbon star counts in great circle cells, llbata et al.l 
(|200lh also found a peak at (90°, 13°) identified with 
Sagittarius. Based on their coui it analyses of M giants 
in the 2MASS survey (A = 5°), iMaiewski et"aLl l}200l 
report poles at (93°, 13°) and (92°, 12°) for the Sagittar- 
ius stream. Dividing their data set into northern and 
southern Galactic hemispheres, they identify poles at 
(99°, 18°) and (91.5°, 11.5°), respectively. 

The choice of width. A, affects the resolution of the 
stream. For the A = 10° adopted in other surveys, the 
highly significant points in Figure [7] coalesce into a sin- 



gle high significance clump. Adopting A < 5° reveals 
additional structure in the stream; however, very small 
A 1 — 2° results in very noisy structures and more 
false positives. If the stream were simple, we could use 
the maximum of the probability pcc as a function of the 
great circle width to determine its approximate width, 
but this approach fails because larger widths pick up 
multiple wraps of Sagittarius. 

Figure E] shows the positions and color-coded radial ve- 
locities for sample stars in the great circles corresponding 
to the 25 poles with the most significant p^fj^^f'^'^ (see 
Table [T]). We also show the position of the Sagittarius 
dwarf galaxy. Most of the highly significant poles are 
associated with the known Sagittarius stream. 

The first set of poles i n Table [1] comprises br anch A of 
the Sagi ttarius stream (Belokurov et aLll2006[ ). In Fig- 
ure 1 of Belokurov et al. ()2006l ). branch A of Sagittarius 
is bounded by a quadrilateral with vertices at approxi- 
mately {a,S) = (125°, 24°), (125°, 17°), (185°, 16°) and 
(185°, 7°). We identify poles corresponding to great cir- 
cles lying entirely within this region with branch A. Our 
two most significant poles thus correspond to branch A. 
The second set of five poles is also associated with Sagit- 
tarius branch A. The great circles for these poles sub- 
stantially overlap the bounding region of branch A. 

A third set of five poles corresponds to other parts of 
the Sagittarius stream. These poles match those found in 
surve ys of carbon stars (llbata et all 1200 1[ ) and M giant 
stars (IMaiewski e"raI1l2003D . Thev_ are also consistent 
with streams identified in 2MASS (|Ibata et al.ll2002D . Al- 
though each of these poles has a smaller p'^^^^'^'^ than 
the most significant poles in branch A of Sagittarius, the 
radial velocity data provide a conclusive detection of halo 
stars associated with the stream. 

Despite these s uccesses, we do not dete ct the Sagittar- 
ius branch B from lBelokurov et al.l (|2006l ) in our 25 most 
significant poles. This failure probably reflects its lower 
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surface density relative to branch A. 

The velocity distributions of stars in the great cir- 
cles corresponding to the poles in the first three sets 
in Table 1 are all similar to the velocity distribution 
shown in Figure [5] with peaks in the velocity interval 
Vrf G [-110, -55] km s"^ 

The last set of five poles in Table 1 includes great cir- 
cles that appear to be superpositions of the Sagittarius 
stream with other structure. Dashed lines in Figure |8] in- 
dicate the corresponding great circles. The first four en- 
tries intersect the Sagittarius stream and the Virgo over- 
density (p^ewberg et al.l 120071 ). Thus, we suspect that 
both of these structures contribute to the overdensities 
of stars in these four great circles. Although the last two 
great circles in Table 1 appear significant in star counts 
and radial velocities, their great circle arc lengths are 
much shorter, 49° and 37°, than those of other great 
circle segments in Table 1. Thus, we cannot robustly as- 
sociate either of these great circles with real structures 
in the halo. 

4. BLUE STRAGGLERS AND BHB STARS IN SAGITTARIUS 

The HVS survey samples a population of blue strag- 
glers and BHB stars that has not been used exten- 
sively to probe the nature of the Sagittarius tidal de- 
bris. From §3.4| the most significant structures in the 
HVS survey result from Sagittarius. We estimate an up- 
per limit on the proportion of stars in the HVS survey 
sample associated with the Sagitarrius stream directly 
from the sample. By comparing the sample of stars in 
o ur most significant strea ms with the A^-body models 
of iLaw fc Majewskl (|201Cll l. we can (1) estimate a lower 
limit on the fraction of stars in the HVS survey sample 
associated with the Sagittarius stream, (2) explore the 
variation in stellar populations along the stream, and (3) 
constrain the surface density of stars within the stream. 

4.1. Stellar Populations and Sagittarius N-body Models 

ILaw fc Maiewskil (j2010( ) construct and analyze a so- 
phisticated A^-body model of the evolution of Sagittarius 
in a triaxial Milky Way potential. The model matches 
all existing observational constraints on the nature of 
the dwarf galaxy and the streams of debris, providing a 
benchmark for interpreting the apparent debris streams 
we detect in the HVS survey. 

We test the positions and radial velocities of HVS sam- 
ple stars in the 25 rn o st sig nificant great circles against 
the ILaw fc Maiewskil (|2Q10D model. Figure M shows the 
position, heliocentric distance, an d galactocentric veloc- 
ity for the lLaw fc Maiewskil ()2010D A^-body model (black 
d ots). 

iBrown et al.l (|2010f ) compute luminosity estimates for 
each star in the HVS survey by matching observed colors 
and spectroscopy to stellar evolution tracks for metal- 
poor mam sequ ence (iGirardi et al.ll2002l [200l and post- 
mam sequence (iDotter et al. 2007, 2008) stars. Coupled 
with the ILaw fc Maiewskil (|201Q) model, these luminosi- 
ties enable us to discriminate between blue stragglers and 
BHB stars in the HVS survey. 

BHB stars in the HVS survey have a mean Mg = 
1.15ib0.35; blue str agglers have a mean Mg = 2.25 ±0.45 
(jBrown et al.ll2010l) . Thus, blue stragglers should sample 
the nearest regions of the Sagittarius stream; BHB stars 
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Fig. 9. — The po sition, h eliocentric distance, and galactocen- 
tric velocity of the ILaw fc M aiewski (2010) Sagittarius A^-body 
model (black dots) and the blue horizontal branch (red squares) 
and blue stragglers (cyan diamonds) that match the Sagittarius 
A'^-body model in position, velocity, and distance. 



should sample the furthest regions. Both types of stars 
contribute to the stream for distances 20 < d < 45 kpc. 

Of the 417 stars contained in the 25 most significant 
great circles, 60 (14%) match the Sagittarius A^-body 
model in position, velocity, and distance for a BHB star 
(red squares in Figure [S]); 51 (12%) match the Sagittar- 
ius A"-body model in position, velocity, and distance for 
a blue stragglers (cyan squares in Figure jS]) ■ Seventeen 
stars are common to both matches. Thus 22.5% of the 
stars in our 25 most significant great circles, roughly 10% 
of the total HVS survey sample, probably belong to the 
Sagittarius stream. Because some stars in the Sagittar- 
ius stream might lie in other great circles, this estimate 
yields a rough lower limit to the fraction of HVS survey 
stars within the Sagittarius stream. 

The HVS survey stars that match the Sagittarius 
stream A^-body model in the region 0° < RA < 50° 
are an equal mix of blue stragglers and BHB stars. The 
stars that match the Sagittarius stream in the region 
130° < RA < 180° are a 2:1 mix of blue stragglers and 
BHB stars (Figure ^ . The HVS survey fairly samples 
both stellar populations within the Sagittarius stream in 
these regions. The observations thus support the idea 
that the stellar population varies along the stream. 

From color-magnitu d e di agrams of Sagittarius, 
iNiederste-Ostholt et al.l (|2010[) find equal numbers of 
blue stragglers and BHB stars in the region 150° < 
RA < 10°. Detailed star count analyses, however, 
suggest a variation in the number ratio of BHB stars 
to main sequen ce turn-off stars along the stream (e.g., 
IBell et all I2010D . Our analysis is also consistent with 
chan ges i n the stellar population along the stream. 
Bell e t al.l (|2010l ) propose that Sagittarius once had 
a BHB-rich core and a BHB-poor halo. Population 
variations then result from stripping different re- 
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Fig. 10. — Estimated number of stars from the HVS survey within 
Sagittarius as a function of great circle width A. 



ffions of the progenitor at di fferent times (se e also , 
Martfnez-Delgado et al. 2 00l IBellazzini et alj 120061: 
Chou et al- ,2007t iPefiarrubia et alJ I2010t ICarlin et atl 
201 in . Because our analysis yields matches to the 
position and velocity of two stellar populations, blue 
stragglers and BHB stars , our resu l ts su pport and 
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extend the conclusion of IBell et al.l ((20101) that the 
stellar population in the Sagittarius dwarf varied with 
radius prior to disruption. 

4.2. Density of Blue Stars in the Sagittarius Stream 

The lLaw fc Maiewskil ()2010[ ) model sets a rough lower 
limit on the contribution of Sagittarius to the HVS sur- 
vey population. From the observed data, we derive an 
approximate upper limit. To make this estimate, we 
compute great circle counts for widths. A, between 1° 
and 20°. For each A, we perform the analysis outlined 
in ij2]and identify the great circle with the smallest p joint- 
Wc then compare the number of stars within this great 
circle, rtofcs, with the number expected, Uexp, for a ran- 
dom distribution of stars on the sky. The estimated num- 
ber of stars within the Sagittarius stream is the excess 
number of stars above random, nsgr = nobs — n^xp- A 
posteriori, we verify that the locations of the poles for 
this collection of great circles change little as a function 
of A and that the excess number of stars as a function 
of A is independent of any small change in the (Z, h) of 
the most significant pole for an adopted A. 

Figure [To] shows our estimate for the number of survey 
stars in the Sagittarius stream, nsgr, as a function of A. 
The derived nsgr rises from ~ 70 for A = 5° to ~ 150 for 
A = 14° and then levels off. Adopting a typical nsgr ~ 
150 for A = 14°-20°, this resuh implies that roughly 
17% (150 stars out of 881) in the HVS survey sample 
belong to the Sagittarius stream. This estimate is larger 
than the ^ 10% derived from the com parison of our 25 
most significant great circles with the iLaw fc Maiewskl 
(pOlO) model in §4.11 Combining these results shows that 
Sagittarius contributes about 10% to 17% of the stars in 
the HVS survey. 

Because this approach uses the full HVS survey, it 
yields a better estimate of the fraction of HVS survey 
stars within the Sagittarius stream. 



4.3. Stellar Surface Density of the Sagittarius Stream 

To conclude this section, we compare the surface den- 
sity of stars identified in our observatio ns of the Sagittar- 
ius str eam with the predictions of the ILaw fc Maiewskl 
(|2010f ) model. For each iteration in a Monte Carlo sim- 
ulation, we construct a sample of model stars with the 
same number of stars observed in the HVS survey. This 
mo del samp le consists of (i) stars drawn randomly from 
the ILaw fc M aie wskil ()2010n model and (ii) stars dis- 
tributed randomly and isotropically in the HVS survey 
area. Based on o ur result in tj4.2i we sele ct 17% of this 
sample from the ILaw fc Maiewskil ()2010[ ) model. We 
then analyze great circle counts using the same procedure 
as for the HVS survey sample with A = 5°. Repeating 
this process 200 times, we derive the median value of 
PGC for the most significant pole in the model Sagittar- 
ius stream. 

Using equation [l] we can express the ratio of the ob- 
served surface density to that predicted from the Monte 
Carlo simulations as 7 = Ih>bsPMc_^ Here, Uohi is the num- 

' riMCPobs ' 

ber of stars in the most significant great circle observed 
in the HVS survey sample, umc is the median value of 
the number of stars in the most significant great circle 
from the Monte Carlo simulations, and Pobs and pmc 
are, respectively, the probabilities that a star placed ran- 
domly in the sample region falls within the observed or 
simulated great circle. 

This analysis shows that the ILaw fc Majewskil ()2010l ) 
model matches our observed surface density for stars in 
the Sagittarius stream. For simulations with A = 5°, we 
derive 7 = 1.03, with a 90% confidence interval of 0.82- 
1.43. Simulations with other A's yield similar results. 

5. CONCLUSION 

We construct a statistical technique to detect star 
streams in the Galactic halo based on both the sur- 
face number density and the radial velocity distribution 
within great circles. The addition of radial velocities to 
the analysis increases our ability to identify star streams 
and to show that some overdensities in number counts 
are superpositions of several distinct structures in veloc- 
ity space. The technique is straightforward yet powerful, 
can be generalized to diverse sets of data, and quantifies 
the statistical significance of tidal star streams. Detec- 
tion is possible even when a structure contributes only a 
few percent of the stars in a sample. 

We apply the method to the HVS survey and detect 
the Sagittarius stream at high significance. Comparison 
with theoretical models confirms that blue stragglers and 
BHB stars are members of the Sagittarius stream. Great 
circle counts and comparisons with theoretical models 
suggest that the Sagittarius stream comprises 10% to 
17% of the halo stars in the HVS sample. The ratio 
of blue stragglers to blue horizontal branch stars varies 
along the length of the stream, with roughly equal num- 
bers of these two stellar types at RA = 0°-50° and a 
2:1 mix of blue stragglers to blue horizontal branch stars 
at RA = 130°-180°. Our conclusions support previous 
indications of variations in stellar population related to 
the original structure of the dwarf galaxy. 

This technique is easily applie d to other large radial 
velocity su rveys including RAVE (iSteinmetz et al.ll2006l: 
IZwitter et al.. .200&) and SEGUE (fVannv et al.l 120091) . 
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Identifying tidal star streams in these radial velocity sur- 
veys will ultimately improve constraints on the Milky 
Way's formation history, dark matter halo mass distribu- 
tion, and the distributions of stellar populations within 
the progenitors. 
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TABLE 1 

The 25 Most Significant Poles Ranked by p^-^.^'^f'"^ 

^ J omt 



7/1 \ 7/1 \ Ti 1 adi II ■■^f ed adjusted 

I (dog) b (clog) Rank par Pv, ^ PJoint P^r 'Pjami 



Sagittarius Stream - Branch A 



101.0 


16.3 


3 


2.29 X 10" 


-12 


1.11 X 


IQ- 


-01 


2.56 X 


10- 


■13 


2.27 


X 


10" 


-08 


2.53 


X 


10- 


-09 


101.6 


15.0 


16 


9.44 X 10" 


-11 


1.62 X 


10" 


-01 


1.53 X 


10- 


■11 


9.33 


X 


10" 


-07 


1.51 


X 


10- 


-07 


102.4 


16.3 


11 


6.61 X 10" 


-11 


4.77 X 


10" 


-02 


3.16 X 


10- 


■12 


6.54 


X 


10" 


-07 


3.12 


X 


10- 


-08 


103.1 


12.4 


1 


3.67 X 10- 


-13 


1.28 X 


10- 


-01 


4.69 X 


10- 


■14 


3.63 


X 


10- 


-09 


4.63 


X 


10- 


-10 


103.7 


16.3 


19 


5.05 X 10- 


-10 


4.92 X 


10- 


-02 


2.49 X 


10- 


■11 


4.99 


X 


10- 


-06 


2.46 


X 


10- 


-07 


103.9 


13.7 


9 


1.35 X 10- 


-11 


1.89 X 


10- 


-01 


2.56 X 


10- 


12 


1.34 


X 


10- 


-07 


2.53 


X 


10- 


-08 


104.3 


15.0 


24 


6.05 X 10- 


-10 


1.17 X 


10- 


-01 


7.08 X 


10- 


■11 


5.98 


X 


10- 


-06 


7.00 


X 


10- 


-07 


104.5 


12.4 


2 


5.05 X 10- 


-13 


2.18 X 


10- 


-01 


1.10 X 


10- 


13 


4.99 


X 


10- 


-09 


1.09 


X 


10- 


-09 


105.3 


13.7 


6 


1.47 X 10- 


-11 


1.12 X 


10- 


-01 


1.65 X 


10- 


12 


1.45 


X 


10- 


-07 


1.63 


X 


10- 


-08 


106.6 


13.7 


15 


4.69 X 10- 


-11 


2.85 X 


10- 


-01 


1.34 X 


10- 


■11 


4.63 


X 


10- 


-07 


1.32 


X 


10- 


-07 



Associated with Sagittarius Stream - Branch A 



96.9 


16.3 


21 


2.55 X IQ- 


-10 


1.19 X IQ- 


-01 


3.04 X 10- 


ll 


2.52 X 10" 


-06 


3.01 X 10" 


-07 


99.6 


16.3 


4 


2.18 X IQ- 


-12 


3.17 X IQ- 


-01 


6.91 X 10- 


■13 


2.15 X 10" 


-08 


6.83 X 10" 


-09 


100.2 


15.0 


25 


7.93 X IQ- 


-10 


9.01 X IQ- 


-02 


7.15 X 10- 


■11 


7.84 X 10" 


-06 


7.07 X 10" 


-07 


101.8 


12.4 


8 


3.13 X 10- 


-11 


7.12 X 10" 


-02 


2.23 X 10- 


■12 


3.09 X 10" 


-07 


2.20 X 10" 


-08 


105.8 


12.4 


10 


1.02 X 10- 


-11 


2.70 X 10- 


-01 


2.76 X lo- 


12 


1.01 X 10" 


-07 


2.73 X 10" 


-08 



Associated with Sagittarius Stream 



82.5 


15.0 


20 


9.98 X IQ- 


-10 


2.61 X 10- 


-02 


2.60 X 10- 


ll 


9.86 X 10" 


-06 


2.57 X 10" 


-07 


83.9 


15.0 


13 


4.20 X 10- 


-10 


1.35 X 10- 


-02 


5.67 X 10- 


-12 


4.15 X 10" 


-06 


5.60 X 10" 


-08 


89.3 


15.0 


12 


2.43 X 10- 


-09 


1.47 X 10- 


-03 


3.56 X 10- 


-12 


2.40 X 10" 


-06 


3.52 X 10" 


-08 


90.7 


15.0 


18 


7.70 X 10- 


-09 


3.08 X 10- 


-03 


2.37 X 10- 


-11 


7.61 X 10" 


-06 


2.34 X 10" 


-07 


91.4 


16.3 


22 


1.99 X 10- 


-08 


2.79 X 10- 


-03 


5.55 X lo- 


-11 


1.96 X 10- 


-04 


5.48 X 10- 


-07 



Superpositions with Sagittarius 



121.6 


20.2 


5 


1.12 X IQ- 


-12 


9.97 X 10" 


-01 


1.11 X 10- 


-12 


1.10 X 10" 


-08 


1.10 X 10" 


-08 


122.1 


18.9 


23 


1.29 X IQ- 


-10 


4.88 X 10" 


-01 


6.31 X 10" 


-11 


1.28 X 10- 


-06 


6.24 X 10" 


-07 


124.9 


30.7 


7 


2.14 X 10" 


-10 


1.01 X 10" 


-02 


2.16 X 10" 


-12 


2.11 X 10" 


-06 


2.14 X 10" 


-08 


126.4 


30.7 


14 


5.21 X 10" 


-10 


1.75 X 10" 


-02 


9.12 X 10" 


-12 


5.15 X 10" 


-06 


9.01 X 10" 


-08 


131.0 


30.7 


17 


1.21 X 10- 


-09 


1.67 X 10" 


-02 


2.01 X 10- 


-11 


1.19 X 10- 


-05 


1.99 X 10" 


-07 



